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We present total and differential cross sections for Wbb and Zbb production at the CERN Large 
Hadron Collider including Next-to-Leading Order (NLO) QCD corrections and full bottom-quark 
mass effects. We discuss the scale uncertainty of the total cross sections due to the residual 
renormalization- and factorization-scale dependence of the truncated perturbative series. We also 
discuss b-quark mass effects in kinematic distributions by comparing with a calculation that con- 
siders massless bottom quarks, as implemented in the Monte Carlo program MCFM. The effects 
of a non-zero bottom-quark mass (m&) cannot be neglected in phase-space regions where the rele- 
vant kinematic observable, such as the transverse momentum of the bottom quarks or the invariant 
mass of the bottom-quark pair, are of the order of m^. Finally, we present the result of a detailed 
comparison of NLO QCD predictions for W + b-jet production with one or two jets with Tevatron 
data. 
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1. Introduction 

The associated production of a weak gauge boson with one or two b jets constitutes not only 
an important background process to Higgs boson searches for light Higgs bosons (M# ~ 135 GeV), 
to single top-quark production and many searches for signals of new physics, but also represents 
a unique opportunity to test and improve the theoretical prediction for heavy quark jets at hadron 
colliders. The cross sections for W and Z boson production with b jets has been measured at the 
Tevatron pp collider at Fermilab (-y/s = 1 .96 TeV) by both the CDF [[T], ||] and DO []|] collaborations, 
and these measurements will continue to improve once more data have been analyzed. Studying 
the same cross sections in the very different kinematic regimes available at the LHC pp collider 
will then be of great interest and will represent a crucial test of our understanding of QCD at 
high-energy colliders. 

The production of a W or a Z boson with up to two jets, one of which is a b jet, has been 
calculated including NLO QCD corrections in the variable-flavor scheme (VFS) or five-flavor- 
number scheme (5FNS) [ffl], while the production of a W or Z boson with two b jets has been 
derived at NLO in QCD using the fixed-flavor scheme (FFS) or four-flavor-number scheme (4FNS), 
first in the massless ft-quark approximation ||] and more recently including full &-quark mass 
effects [0, j^,^, 10]. In the FFS or 4FNS only massless-quark densities are considered in the initial 



state, while in a VFS or 5FNS an initial-state ft-quark density is introduced. The two schemes 
amount to a different ordering of the perturbative series for the production cross section: in the 
4FNS the perturbative series is ordered strictly by powers of the strong coupling a s , whereas in the 
5FNS the introduction of a b-quaik partem distribution function (PDF) allows to resum terms of 
the form a"\n{m 2 h /M 2 ) m at all orders in a s (for fixed order of logarithms m), where M represents 
the upper integration limit of the &-quark transverse momentum and can be thought to be of the 
order of Mw or Mz- While the two approaches can give very different results at leading order (LO) 
in QCD, starting at NLO in QCD the total cross sections have been shown to be consistent within 



their respective theoretical uncertainties for both H+l ft-jet production [11] (for a brief review see 



also Ref. [12]) and single-top production [|l3|]. In kinematical distributions, however, &-quark mass 
effects which are fully taken into account in the 4FNS can have a significant impact, particular in 
phase space regions where the relevant kinematic observable is of the order of m\,. Bottom-quark 
mass effects are also important in W + 1 £»-jet production and improved predictions at NLO QCD 



are provided in [14], where the calculations of Wbb production in the 4FNS with massive b quarks 



and Wbj production in the 5FNS have been merged. A similar study is currently in progress for 



Z + 1 ft-jet production [15]. Improving the predictions for Z + 1 b-jet production will be particularly 
relevant at the LHC, where this process allows for a direct determination of the b-quaik PDF, to 
be used in the prediction of H + 1 ft-jet production, a discovery channel for beyond-the-SM Higgs 
bosons with enhanced ft-quark Yukawa couplings. 

In these proceedings, we present results for Wbb and Zbb production at the LHC, i.e. with 
both b jets tagged in the final state, keeping the W and Z boson on shell. We discuss the impact of 
the ft-quark mass on total kinematic distributions by comparing with a calculation which neglects 
the fe-quark mass at NLO in QCD, as implemented in the MCFM package JT6[] . Finally, we present 
the results of a detailed comparison of NLO QCD predictions for W + &-jet production with one or 
two jets with Tevatron data. 
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2. Wbb and Zbb production at NLO QCD 

At tree level, the production of a W boson with a pair of bottom quarks consists of just one 
process, qq' — > Wbb, while Zbb production consists of two channels, namely qq — > Zbb and gg — > 
Zbb. In order to compute these processes at NLO in QCD one needs to include one-loop virtual 
corrections to the tree-level processes, as well as all real radiation corrections with up to one extra 
partem in the final state, i. e. qq' — > Wbb + g, gg,qq —> Zbb + g, and the ^-initiated processes, 
qgiqg) ~^ Wbb + q' {q') and qg(qg) —> Zbb + q{q). Details of the calculations and numerical results 
are provided in Refs. [@, |, g, go||. 

The LO results are based on the one-loop evolution of a s and the CTEQ6L1 set of PDFs [|l7|], 
with af°(M z ) = 0.130, while the NLO results use the two-loop evolution of a s and the CTEQ6M 
set of PDFs, with a^ L0 {Mz) = 0.118. We use the kj jet algorithm with a pseudo-cone size of 
R = 0.7 and we checked that our implementation of the kj jet algorithm coincides with the one in 
MCFM. We require all events to have a Z?-jet pair in the final state, with a transverse momentum 
either larger than 15 GeV or 25 GeV and require the pseudorapidities of both Z? jets satisfy \r\ b ' h \ < 
2.5. We impose the same pj and |t]| cuts also on the extra jet that may arise due to hard non- 
collinear real emission of a parton, i. e. in the processes W /Zbb + g or W /Zbb + q{q). This 
hard non-collinear extra parton is treated either inclusively or exclusively. In the inclusive case we 
include both two- and three-jet events, while in the exclusive case we require exactly two jets in the 
event. Two-jet events consist of a &-jet pair that may also include a final-state light parton (gluon 
or quark) due to the applied recombination procedure. On the other hand, three-jet events consist 
of events containing a ft-jet pair plus an extra light jet. We notice that, at NLO in QCD, all jets in 
three-jet events consist of a single parton. 

In Fig. (W ± bb) and Fig. | (Zbb) we illustrate the renormalization- and factorization-scale 
dependence of the LO and NLO total cross sections obtained for a massive b quark, when /I = 
\l r = \if is varied between Ho/2 and 2/j.q, with /lo = My + 2nib (V = W,Z). We immediately 
notice that the impact of NLO QCD corrections is very large, in particular in the inclusive case 
in Wbb production, where they increase the LO cross section by a factor between two and three 
depending on the scale. We also notice that the scale dependence of the NLO Wbb cross section 
is worse than (inclusive case) or comparable to (exclusive case) the scale dependence of the LO 
cross section. This is different from what has been observed for the Tevatron [f7|, |9|], and was first 
pointed out in a calculation with massless bottom quarks [|6|] . It is j ust a reminder of the fact that, at a 
given perturbative order, the uncertainty due to the residual renormalization- and factorization-scale 
dependence may underestimate the theoretical uncertainty due to missing higher-order corrections. 
A realistic determination of this uncertainty is usually much more complex and requires a thorough 
understanding of the perturbative structure of the cross section, in particular at the lowest orders of 
the perturbative expansion. In both Wbb and Zbb production the NLO QCD corrections introduce a 
new production channel not present at LO, the ^-initiated processes, and as such introduces a LO 
scale dependence. The NLO Wbb total cross section is particularly affected by this process because 
there is no gg-initiated process at LO. In Zbb production, however, this process is not dominant at 
NLO and its impact is therefore less pronounced. Indeed, the scale dependence of the exclusive 
Zbb cross section actually greatly improves at NLO in QCD, while the scale dependence of the 
inclusive one is only mildly better than at LO, but not worse as it is the case in Wbb production. 
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Figure 1: Dependence of the LO (black, solid), NLO exclusive (blue, solid), and NLO inclusive (red, solid) 
total cross sections for W + bb and W~bb production on the renormalization/factorization scales at the LHC 
(^/s = 14 TeV), including full b-quark mass effects, when ji = /!, = /!/ is varied between Hq/2 and 2/Xq 
(with /Xo = Mw + 2nib). We also show the individual parton-level channels, qq' (dash-dotted) and qg + qg 
(dotted), for the inclusive (red) and exclusive (blue) cases. Taken from Ref. [nq]. 



Finally, the impact of the ^-initiated processes on the scale dependence of the total cross sections 
is larger in the inclusive than in the exclusive case because the exclusive cross section by definition 
discriminates against processes with more than two jets in the final state. 

In Fig. [| we illustrate ft-quark mass effects on the example of m h i distributions in W~bb and 
Zbb production. These effects can impact the shape of the kinematic distributions in particular in 
phase space regions where the relevant kinematic observable is of the order of nit,. These effects can 
be approximated by rescaling the NLO cross section for nib = with the ratio of LO cross sections 
for massive and massless bottom quarks as discussed in detail in [J7| @]. The total production cross 
sections are reduced by ft-quark mass effects, and the effect is more pronounced the smaller the 
applied p b T cut. However, these effects are in most cases smaller than the residual scale dependence 
at NLO in QCD, especially in Wbb production for the inclusive case. 

3. Comparison of W + b-jet predictions with Tevatron data 

A recent measurement of the cross section of W boson production in association with one or 
two b jets by the CDF collaboration at the Tevatron finds [Q] 

Ofc-jeu x 3S(W -»■ ^v)(CDF) = 2.74±0.27(stat.) ±0.42(syst.) pb 

This ft-jet cross section includes W + 2£-jets and W + lZ?-jet contributions, and the NLO QCD 
predictions for both signatures have to be considered. The NLO QCD prediction for W + 2^-jets 
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Figure 2: Dependence of the LO (black, solid), NLO exclusive (blue, solid), and NLO inclusive (red, 
solid) total cross sections for Zbb production on the renormalization/factorization scales at the LHC (y/s = 
14 TeV), including full fo-quark mass effects, when /x = \i r = fif is varied between jUo/2 and 2/io (with 
jUo = Mz + 2m/,). We also show the individual parton-level channels, qq' (dash-dotted), qg + qg (dotted) and 
gg (dashed), for LO and NLO (inclusive (red) and exclusive (blue)). Taken from Ref. Jl(i|]. 



production is based on Refs. [{7|, |9|] and the one for W + lb-jet production on Ref. [jfj], where in 
both cases events with a non-Z? jet that result in a three-jet event are discarded. In Ref. [|l4|], the 
calculations of Wbb production in the 4FNS and Wbj production in the 5FNS have been merged 
by consistently treating contributions which arise in both processes. In this way, the Z?-quark mass 
is taken into account, which cannot be neglected when one b quark is treated fully inclusively, the 
tree-level ^-initiated process discussed earlier is improved by including NLO QCD corrections 
in the 5FNS, and large initial-state logarithms are resummed via the b-quaik PDF approach. As a 
result of this improved treatment, theoretically more stable predictions for W + lZ?-jet production 
are now available. Finally, combining the results of predictions for W + 2Z?-jets and W + lft-jet 
production, thereby doubling the weights of 2^-jet events to obtain jet-level cross sections, yields 



the following NLO QCD prediction (with \i r = [if = M w for the central value) []18p: 



cruets x SS{W -> ^v)(NLO QCD) = 1.22±0.14(scale uncert.) pb . 

Together with the LO prediction of 0.91^0 20 (including scale uncertainties) this results in a 
moderate K-factor of about 1.34. The origin of the discrepancy between theory and experiment is 
presently under investigation 1 . 



Acknowledgments 

We would like to thank the organizers of the workshop Northwest Terascale Research Projects 
W + b quark physics at the LHC, especially Davison Soper, for the kind hospitality extended to 

1 See also talks given at the workshop Northwest Terascale Research Projects W + b quark physics at the LHC, held 
at the University of Oregon, http://physics.uoregon.edu/-~soper/TeraWWW2. 



5 



Associated production of a W or Z boson with bottom quarks 



Doreen Wackeroth 



20 40 60 80 100 120 140 160 ISO 200 Q 20 40 60 80 100 120 140 160 180 200 

1 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 10 C 1 I ' I ' I ' | ' I ' | ' I ' | ' I 1 I 




Figure 3: LO (black), NLO inclusive (red) and NLO exclusive (blue) m hb distributions for Zbb (left hand 
side) and Wbb (right hand side) production at the LHC derived with m& 7^ (LO: dashed, NLO inclusive: 
solid, NLO exclusive: dash-dotted) and from MCFM with mt, = (LO: double-dashed/dotted, NLO /nc/w- 
sive: dashed/double-dotted, NLO exclusive: dotted). The lower window shows the ratio of the distributions 
for massive and massless b quarks, da[m\, 7^ 0)/da(mi, = 0) (LO: dashed, NLO inclusive: solid, NLO 
exclusive: dash-dotted). Taken from Ref. jio|]. 
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